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Abstract 

Recently it was shown that the eigenfunctions for the the asym- 
metric exclusion problem and several of its generalizations as well 
as a huge family of quantum chains, like the anisotropic Heisenberg 
model, Fateev-Zamolodchikov model, Izergin-Korepin model, Suther- 
land model, t — J model, Hubbard model, etc, can be expressed 
by a matrix product ansatz. Differently from the coordinate Bethe 
ansatz, where the eigenvalues and eigenvectors are plane wave combi- 
nations, in this ansatz the components of the eigenfunctions are ob- 
tained through the algebraic properties of properly defined matrices. 
In this work, we introduce a formulation of a matrix product ansatz for 
the six-vertex model with periodic boundary condition, which is the 
paradigmatic example of integr ability in two dimensions. Remarkably, 
our studies of the six- vertex model are in agreement with the conjec- 
ture that all models exactly solved by the Bethe ansatz can also be 
solved by an appropriated matrix product ansatz. 

1 Introduction 

The Bethe ansatz in its several formulations (coordinate, inverse scattering 
and functional) has been established over the years as a powerful tool for the 
description of the eigenvectors of a huge variety of integrable one-dimensional 
quantum spin chains and two-dimensional transfer matrices (see e.g. PQ-[S] 
for reviews). A quantum Hamiltonian or a transfer matrix is considered ex- 
actly integrable if an infinite number of its eigenstates can be expressed by 
the Bethe ansatz in the thermodynamic limit. In the last two decades it has 
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been shown that a matrix product ansatz (MPA) can be used to express the 
stationary distribution of the probability densities of some special stochastic 
models [6]- [9]. Although these models are in general not integrable through 
the Bethe ansatz, they have the components of its ground-state wavefunc- 
tions given in terms of a product of matrices. According to this ansatz, the 
algebraic properties of the matrices defining the MPA fix the components of 
the wavefunction apart from a normalization constant. 

An important development of the MPA that appeared in the context of 
stochastic models is the dynamical matrix product ansatz (DMPA) [TU1 [TT] . 
This DMPA was shown originally to be valid to the problem of asymmetric 
diffusion of particles on the lattice |10j and extended to other stochastic mod- 
els and related spin Hamiltonians [T2| [T3] . This ansatz allows the calculation 
of the probability densities, of the stochastic system, at arbitrary times. In 
the related spin Hamiltonian this DMPA asserts that not only the ground- 
state wave function, as in the standard MPA, but an arbitrary wavefunction 
have its components expressed in terms of a matrix product ansatz whose 
matrices, in distinction of the standard MPA, are now time dependent. 

More recently [T3J-JTB] it was shown that several exactly solvable Hamil- 
tonians, related or not to stochastic models, may also be solvable by an 
appropriate time independent matrix product ansatz. In this new MPA not 
only the ground-state but all wavefunctions can be expressed by a product 
of matrices. Using this new MPA it was possible to rederivethe results pre- 
viously obtained through the Bethe ansatz for several quantum chains with 
one and two global conservation laws, such as the XXZ chain, spin-1 Fateev- 
Zamolodchikov model, Izergin-Korepin model, Sutherland model, t-J model, 
Hubbard model, etc (THUS], as well as the exact solution of the asymmetric 
exclusion problem with particles of arbitrary size [TU]. Moreover, the com- 
ponents of the eigenfunctions of the exact integrable Hamiltonians, which 
according to the Bethe ansatz are normally given by a combination of plane 
waves, can also be obtained from the algebraic properties of the matrices 
defining the new MPA. In the case of Bethe ansatz solutions the eigenvalues 
and the amplitudes of the plane waves are fixed apart from a normalization 
constant by the eigenvalue equation of the Hamiltonian. On the other hand, 
in the new MPA, the eigenvalue equation fix the commutation relations of the 
matrices defining the ansatz. The advantage of the new MPA in the search 
for new exact integrable models, as showed in our previous works [Il]-[l6j, is 
its simplicity and unifying character in the implementation for arbitrary sys- 
tems. All the previous successful applications of the new MPA [E]-[IS! was 
concerned with the eigenspectrum of quantum Hamiltonians and stochas- 
tic models. In this paper we are going to show that these results can also 
be extended to transfer matrix calculations of two-dimensional classical spin 
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models. More specifically we are going to extend the new MPA introduced in 
[E]-[IS] for the case of the row-to-row transfer matrix of the six- vertex model 
with toroidal boundary condition, which is the paradigm of integrability in 
two dimensions. This transfer matrix was diagonalized through the coordi- 
nate Bethe ansatz firstly by Lieb [U], in a special case, and by Sutherland 
[IH] and Yang [TIED!. 

2 The asymmetric six-vertex model and its 
transfer matrix 

The six-vertex model defined on a square lattice, was introduced to explain 
the residual entropy of the ice [I7j-[20j. We are going to consider the asym- 
metric version of the six-vertex model that was first studied and exactly 
solved with standard methods [2T1 [22] . This model is defined on a square 
lattice with M rows and N columns and toroidal boundary condition. At 
each horizontal (vertical) lattice bond we attach an arrow pointing to the 
left or right (up or down) directions. These arrows configurations can be 
equivalently described by the vertex configurations of the lattice. A vertex 
configuration at a given site (center) is formed by the four arrows attached 
to its links. The allowed vertex configurations are those satisfying the ice 
rules: two of the arrows pointing inward and the other two pointing outward 
of its center. There are six possible configurations for the vertices. Theses 
configurations are showed in fig. la. In fig. lb, a more convenient notation 
is introduced, in which we only draw by a solid line (broken line) the links 
having arrows pointing to the left or down (right or up) of the center defin- 
ing the vertex. Labeling the M rows sequentially by m = 1, 2, M and by 
{x m } the solid line positions on the vertical edges of the row m, the partition 
function can be written as 

Z = E E • • • E {* 2 }) T ({* 2 }> {* 3 }) ■ ■ • T ({x M }> {z 1 }) = Tr(T M ), 

KHz 2 } {x M } 

(1) 

where T is the 2 N x 2 N transfer matrix, with elements 

T({y}, {x}) = e" /3(niei+n2e2+ - +n6£6) , (2) 

where the summation is over all allowed arrangements of lines on the hor- 
izontal edges and rij (j = 1,...,6) are the numbers of vertices of types 
(1, ...,6) formed by the configurations. For convenience we label the Boltz- 
mann weights ao, ai, W, b 2 , ci, associated with the vertices as in fig. 1. 
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It is also important to mention that the number of vertical and horizontal 
lines is conserved, forming continuous non-crossing paths through the lat- 
tice. On the other hand the transfer matrix, due to the toroidal boundary 
condition, is translation invariant. As a consequence of these symmetries 
the transfer matrix breaks up into blocks of disjoint sectors labeled by the 
number n of vertical lines (n = 0, N) and the momentum eigenvalues P 

(p = fU = o,i,...,iv-i). 

3 The Matrix Product ansatz for the six- vertex 
model 

The ansatz we propose [H]-[16] states that any eigenfunction \^ n ,p) of the 
transfer matrix ([2]) in the sector with n (n — 0, 1, 2, . . . , TV) vertical lines and 
momentum P (P — = 0,1, ...,N — 1) is given in terms of a matrix 

product, i. e., their amplitudes are given by the trace of the following matrix 
product: 

* 

|Vvp>= Yl ^(E xl ~ l AE X2 - xl - 1 A---E* n - Xn - 1 - 1 AE L - Xn n P )\x 1 ,...,x n ), 

(3) 

where \x\, . . . ,x n ) denote the configurations with vertical lines at positions 
(xi, . . . , x n ) and the symbol (*) in the sum means the restriction to the con- 
figurations where L > > x^ > 1. The objects A, E and Qp are abstract 
matrices, or operators, with an associative product whose commutation re- 
lations will be fixed by imposing the validity of the eigenvalue equation of 
the transfer matrix (j2J). The matrices A and E are associated with the sites 
where we have a vertical line or not, respectively, and the matrix Qp is in- 
troduced in order to fix the momentum P of the eigenfunction |\l/ n p). The 
fact that \ipn,p) has a momentum P imply that the ratio of the amplitudes 
corresponding to the configurations \x\, . . . , x n ) and \x\ + 1, . . . , x n + 1) is 

— iP ■ 

e , i.e., 



Tr( E^- 1 ■AE* 3- * 1-1 A ■ ■ ■ E Xn ~ x "~ 1 ~ 1 AE L ~ Xn Q P ) 



(4) 



Tr(E x ^AE x »- x ^- 1 A ■ ■ ■ E x ™- x ^- l AE L - x ™- l Q, P ) 

and consequently from ([3]) we obtain the following commutation relations 

AVLp = e~ iP n P A, EVLp = e- iP tt P E. (5) 

The matrix product ansatz will be valid if the algebraic relations among 
the matrices A, E and Qp are consistent with the constrains imposed by the 
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eigenvalue equation 

T\ip n ,p) = K\^n,p)- (6) 

To solve this last equation, it is helpful to begin, as usual, by considering the 
simple cases where n — 0, 1 and 2 before considering the general case. 
The case n = 0. 

In this case the solution of the eigenvalue equation © is trivial since we 
do not have vertical lines between two successive rows. There are only two 
possible horizontal arrangements either all bonds have a line or all of them 
are empty. In this case the vertices are all of type 1 or type 3 (see figure 1) 
and consequently the eigenvalue is given by 

Ao = 4 + &f (7) 

where ao and b\ are the Boltzmann weights of the vertices of types 1 and 3, 
respectively (see figure 1). 
The case n = 1. 

We have in this case just one vertical line between two rows. The transfer 
matrix links a vertical line at position y (y = 1, ..,7V) above a row to a 
vertical line at any position x (x = 1, ...,7V) under this row. The elements 
of the transfer matrix T(y,x) in this sector with momentum P are given by 
([2]) • They are the product of the Boltzmann weights of the vertex appearing 
on the row. If the position of the line x is less (greater) than y, the vertex 
configuration at these sites will be of types 5 and 6 (6 and 5) and all the 
others vertices will be of types 3 (1) and 1 (3) depending on whether the 
vertices are between the positions x and y, or not, respectively. In the case 
where x = y these vertices should be of type 4 or 2 with all the remains 
vertices of type 1 or 3 respectively. Consequently the eigenvalue equation 
for the transfer matrix ([2]) associated with the components of \ip n ,p) © with 
n = 1 and momentum P give us the relations 

N 

X{lT(E x - 1 AE N - x tt P ) = a^- y+x ' 1 b v 1 - x - 1 c l c 2 i:r(Ey- 1 AE N -yQp) + 

y=x+l 
x-1 

al- v - l b^- x+y - l c 1 c 2 Tr{E y - 1 AE N - y n P ) + 

y=l 

(a$-% + 6f - l ai)Tr (E x ~ l AE N ~ x Qp) . (8) 

Equation ([8]) can be simplified in order to express all the matrix products in 
terms of a single one. This is possible by exploring the cyclic property of the 
trace as well as the commutation relations This allow us to factorize the 
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matrix product in the following form 

N x-1 
y=x+l 3/=l 

(a%- 1 b 2 + b?- 1 a 1 ). 
By evaluating the sums in © we obtain 

r r f h \ ~ X p -iP(l-x) 



^ = a»L{P) + b»M{P) + br—[-\ , iP (l-e-^), (10) 



where 



«o&2 + (cic 2 - &i&2)e lP , , a§a\ - c\c 2 - a x b\e %r 



L(P)= u z ' v \ V and M(P)_ , , , „ 

(11) 

In order to satisfy (jSJ), the eigenvalue Ai in ffTUj) should be independent of 
the vertical line position x. Thus the last term in the right hand side of 
( jTUl) must vanish. The only way to cancel this term, for non zero Boltzmann 
weights, is obtained by imposing the following constraint to the momentum 
P 

e lNP = 1, (12) 

which is automatically satisfied, since P = j^-l I = 0,1,..., AT — 1. The 
eigenvalue (JTUl) is then given by 

A 1 = agL(P)+b?M(P). (13) 

An alternative solution of ([H]), whose generalization will be convenient for 
arbitrary values of n, is obtained by expressing the matrix A in terms of the 
matrix E and a spectral parameter dependent matrix 

A = A k E, (14) 

with A k satisfying 

EA k = e ik A k E. (15) 
As a consequence of © and (jT5"|) A k also satisfies 

A k VL P = VL P A k . (16) 

The spectral parameter k will be fixed by the eigenvalue equation (151) . In- 
serting (j!4p in (jSJ) and using the commutation relation 0151) we obtain © 
with the value k replacing P. Therefore 

Ai = a$L(k)+b?M(k), (17) 
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with 

e w * = l, k = ^l (Z = 0,l,...,iV-l). (18) 

Comparing ( T121) and ( lT3i) with (JTTj) and ( |T8l) we observe the equality k = P. 
This fact can also be seen directly by inserting ( fl4l) in (Tj0) and using ( 1151) . 

We still need to verify whether the algebraic relations among the matrices 
A k , E and Qp (TT3T) and (|TB|) are consistent with the cyclic property of 
the trace. Indeed these equations yield 

Tr(A k E N Q P ) = e- iNk Ti(E N A k n P ) = e~ iNk Ti(E N Q P A k ) 

= e- iNk Tr(A k E N n P ), (19) 

which satisfies the cyclicity of the trace due to ffTSl) . Since no new constraints 
is obtained for the matrices A k , E and Qp, with k = P, and for spectral 
parameter k, the MPA is consistent. 
The case n = 2. 

In this sector there are two vertical lines in the row. We have in general 
two types of relations, which are relations where at least one of the vertical 
lines (1/1,1/2) coincide with (xi,x 2 ) and those where y\ and y 2 interlace with 
X\ and x 2 (xi < Hi < x 2 < IJ2 or y\ < x\ < yi < X2). Then, the eigenvalue 
equation (J6J) imply 

A 2 Tr(E x '- 1 AE X2 - Xl - 1 AE N - X2 Q P ) = 

X2 N * 

E E ^- y ^-\ 2 f{x liyi )g{y u x 2 )f{x 2 ,y 2 )Tr{Ey^ + 

2/1=2:1 V2=X2 
X\ X2 * 

E E K~ x ^- l c ig { yi , x 1 )f(x 1 , y 2 )g(y 2 , x 2 )Tr(E yi ~ 1 AE y2 ~ yi ~ 1 AE N ~ y2 Q P ){20) 

?/l=l 2/2=3:1 

where the symbol * in the sums means that terms with yi = y 2 are excluded 
and 

*, s f - if a; = 2/ , , N f ^ if x = y 

y) = 1 cxftT"- 1 if » > * ^ ^ (?/ ' X) = I c 2 af - 1 if x > y • 

(21) 

The relation (1201) connects configurations where the arrangements of vertical 
lines above one row do not have the same distance of the vertical lines below 
this same row. In other words, the distance of the incoming lines y 2 — y\ are 
in general different of the outcoming distance x 2 — x\. As a consequence, 
it is not possible to solve the eigenvalue equation by just using the cyclic 
property of the trace in as done previously in the case n = 1. We need 
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now to use a generalization of the algebraic relation (|14|) for the case of two 
lines. The generalization of (|14p is done by writing the matrix A in terms of 
two new spectral parameter matrices as: 



3=1 



'.22) 



with the commutation relations 

EA k] =e ik ^A kj E and A kj Q P = Q P A kj (j = 1,2), (23) 

where the spectral parameters k\ and k 2 are up to now unknown complex 
numbers. 

Inserting fl22|) in (1201) and using in this expression (123]) and ffTTj) we ob- 
tain, after similar manipulation as we did in the case n — 1, the following 
constraints 



[A 2 - afiUkJUki) - b? MfaWfa)] e- lk ^e- lk ^Ti(A k A kl E N n P ) 

3,1=1 

2 



3,1=1 

-£»f 

3,1=1 

3,1=1 



L{ki)M{kj) 



L(^)M(%) - 



a x b 2 




a bi_ 




aib 2 




a b 1 _ 


(oo 



-i(kj+ki)x 2 



Tr(A k .A kl E"n P ) 



X\—X2 



e-^+^Tr (A k .A kl E N Q P ) 



C 2 C 2 |~g— iNk tp — ikjX 



g lKjX± g 



-ikix 



i^ij e -i{kj+ki) / ^ 



Qo( a o — b\e tkj )(a — b\e~ %k 



-J- 2 



-± Ti(A kj A kl E»n P ) 
a 0- 



i,i=i 

2 



C 2 C 2 fg— t(iV+l)fcj g — ifcj-ara _ g-iAy'g-ifei^J / ^ 



+ ^ 6f cic 2 6 2 

3,1=1 
2 

b x c\c 2 a\ 

3,1=1 



a bi(a — 61 e ifej )( a o — &ie~* fc ') 

— i(iV+l)fc; g—ikjXi 



- Tr(A kj A h E»n P ) 



e J 6 



*o( a o - 6ie~ ifc; ) ag(a -6ie- ifc 



a 



-i(N +l)ki p —ikjX2 



iki „—ik,X2 



-e ] e 



a bi(a — b 1 e~ lkl ) a bi(a — b\e %k i~) 



Tr(A k] A kl E"n P 
h x - Xl 



a 



(24) 



Tr^.A^^Op) = 0, 



where 1 < X\ < x 2 < A/". This can only be satisfied if each sum is identically 
zero. Moreover since A 2 should be independent of x\ or x 2 a possible solution 
of (Ell) i s obtained by imposing 



A 2 = a%L{k x )L{k 2 ) + b l lM{k x )M{k 



N ; 



(25) 
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The algebraic relation between the matrices A kl and A k2 are obtained by 
imposing that both the second and third terms in (l2~ij) are zero independently, 



A kj A kl = -S{k h h)A kl A k . (l^j) (A k .) =0 (j, Z = 1,2), (26) 
where 



L(k)M(k j; aobi 



s(k j ,h)= - ^ , (27) 



with and M(k) given by (II ip . Finally, the vanishing of the last four 
terms in ( 1!M1) will give us relations that fix the spectral parameters values 
£4 and k 2 . These equations are obtained by exploring the algebraic relations 
(J26D 

e ^h = _ S{ k. M ) (lj = l }2 and l^j). (28) 

The eigenvalues and eigenvectors are obtained by inserting the solutions 
(ki, k 2 ) of these last equations in ([2"o"]) and (12"7|) . respectively. The momentum 
P is obtained by using $22]) and ([23]) in (jlj), i.e., P = h + k 2 . 

The consistency of the algebraic equations ( 1221) . f l23l) and ( 1271) with the 
cyclic property of the trace in ([3]), as in the case n = 1, can be easily verified, 
yielding 

Tr(A k .A kl E N n P ) = -Sfa, h)Tr(A kl A k .E N n P ) 

= -S(k v h)e- lNk ^i(A h E N A k n P ) 

= -S(kj, h)e- iN ^Tr(A k .A kl E N n P ). (29) 

The case of general n. 

The previous calculation can be extended for arbitrary values of the num- 
ber n of vertical lines. The transfer matrix ([2]) when applied to the amplitudes 
of \ip n ,p) give us an eigenvalue equation linking an arrangement of vertical 
lines xi, ...,x n with arrangements yi,...,y n with x\ < y\ < x 2 < • • • x n < y n 
and yi < x% < y 2 < • ■ ■ y n < x n . To solve this eigenvalue equation we need 
to extend the definition ( 1221) and the commutation relations ( 1231) for general 
n, i. e., 

n 

A = J2 A k,E (30) 

with 

EA k . = e ik >A k] E and A k] Vt P = Vt P A k] (j = 1, n), (31) 

where kj (j = 1, . . . , n) are in general unknown complex numbers that will 
be fixed by the eigenvalue equation ([6]). Actually the definition (1301) is not 
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the only possible one. The most general definition that enable us to solve 
the eigenvalue equation is A = Y^ij=i E a A kj E 1+l3 , where a and (3 are integer 
numbers. However ( !30l) is more convenient since otherwise the S-matrix in 
( |32l) and the Bethe equation ( 1331) will depend on the parameters a and f3. 
Inserting ( 1301) in the eigenvalue equation (jBJ) and using the commutation 
relations (1311) we obtain, similarly as done in the case n = 2, the algebraic 
relations among the matrices {Afc.} 

A k .A h = -S(k 3 ,k)A h A k] (A k f = {j±l = l,...,n), (32) 

where S(kj, ki) is given by ( 1271) and the spectral parameters kj (j = 1, ...,n) 
are fixed by the equation 

n 

e iN h = { _ lT -x "Q s{k . )k) ( i = i 7 ... ;n) . (33) 
J=i (#i) 

No new algebraic relations appear for the matrices {A^ . } and the associativity 
of the algebra (|3T|) and (}3"2"1) follows from the property S(kj, ki)S{k h kj) = 1. 
The eigenvalues for the transfer matrix in the sector with general n is 
then given by 

A 2 = a^L^Lih) ■ ■ ■ L{k n ) + b^M{k x )M{k 2 ) ■ ■ ■ M(k n ), (34) 

where L(k) and M(k) are given by (jTTj) and the spectral parameters {kj} 
are the solutions of (1HH1) . The eigenvalues (IMI) and the spectral parameter 
equations coincide with the corresponding equations obtained through the 
Bethe ansatz 

Finally, the momentum P follows from (j4]) and fl3~T|) : 

n 

The consistency of the algebraic relations of the matrices defining the 
MPA OH]) with the cyclic property of the trace in is promptly verified as 
in the cases where n = 1 and n = 2. Therefore the MPA is consistent and a 
infinite number of eigenvectors of the transfer matrix (T5]) can be written by 
([3]) in the thermodynamic limit. 

4 Conclusion 

In conclusion, we have shown that the new MPA introduced in [HJ [T5] for 
one dimensional quantum spin chains, such as the XXZ chain, spin-1 Fateev- 
Zamolodchikov model, Izergin-Korepin model, Sutherland model, t-J model, 
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Hubbard model, as well as the exact solution of the asymmetric exclusion 
problem with particles of arbitrary size [16J, can also be extended to the 
diagonalization of the row-to-row transfer matrix of the six-vertex model 
with toroidal boundary condition. Differently from the standart MPA [6]- [9] 
this new MPA [T4"| [T5] asserts that all wavefunctions can be expressed by a 
product of matrices. The solution of the six vertex model through the new 
MPA is in agreement with the conjecture proposed in [TU] and [THUS] that all 
models exactly solved by Bethe ansatz can also be solved by an appropriate 
MPA. An interesting problem for the future is the formulation of a MPA for 
others spin models like the 8- vertex model, which is related to quantum spin 
chains with no global conservation laws such as the XYZ chain. 
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S. Sarandy for reading the manuscript. This work has been supported by 
CAPES and FAPESP (Brazilian agencies). 

References 

[1] R. J. Baxter, Exactly solved models in statistical mechanics (Academic 
Press, New York, 1982). 

[2] G. M. Schiitz, Integrable Stochastic Many-body Systems in "Phase 
Transition and Critical Phenomena", vol. 19, Eds. C. Domb and J. 
L. Lebowitz (Academic, London, 2000). 

[3] V. E. Korepin, A. G .Izergin and N. M. Bogoliubov, Quantum Inverse 
Scattering Method, Correlation Functions and Algebraic Bethe Ansatz 
(Cambridge University Press, Cambridge, 1992). 

[4] F. H. L. Essler and V. E. Korepin, Exactly Solvable Models of Strongly 
Correlated Electrons (World Scientific, Singapore, 1994). 

[5] P. Schlottmann, Int. J. Mod. Physics B, 11, 355 (1997). 

[6] B. Derrida, M. R. Evans, V. Hakim and V. Pasquier, J. Phys. A: Math. 
Gen. 26, 1493 (1993). 

[7] M. R. Evans, D. P. Foster, C. Godreche and D. Mukamel, J. Stat. Phys. 
80, 69 (1995). 

[8] B. Derrida, J. L. Lebowitz and E. R. Speer, J. Stat. Phys. 89, 135 (1997). 



11 



[9] F. C. Alcaraz, S. Dasmahapatra and V. Rittenberg, J. Phys. A: Math. 
Gen. (N.Y.) 31, 845 (1998). 

[10] R. B. Stinchcombe and G. M. Schlitz, Phys. Rev. Lett. 75, 140 (1995); 
Europhys. Lettt. 29, 663 (1995). 

[11] G. M. Schiitz, Phase Transitions and Critical Phenomema vol. 19 ed C. 
Domb and J. Lebowitz (London: Academic, 2000). 

[12] T. Sasamoto and M. Wadati, J. Phys. Soc. Japan 66, 2618 (1999). 

[13] V. Popkov, M. E. Fouladvand and G. M. Schiitz, J. Phys. A 35, 7187 
(2002). 

[14] F. C. Alcaraz and M. J. Lazo, J. Phys. A: Math. Gen. 37, LI (2004). 

[15] F. C. Alcaraz and M. J. Lazo, J. Phys. A: Math. Gen. 37, 4149 (2004). 

[16] F. C. Alcaraz and M. J. Lazo, Braz. J. Phys. 33, 533 (2003). 

[17] E. H. Lieb, Phys. Rev. 162, 162 (1967). 

[18] B. Sutherland, Phys. Rev. Lett. 19, 103 (1967). 

[19] C. P. Yang, Phys. Rev. Lett. 19, 586 (1967). 

[20] C. N. Yang, Phys. Rev. Lett. 19, 1312 (1968). 

[21] I. M. Nolden, J. Stat. Phys. 67, 155 (1992). 

[22] D. J. Bukman and J. D. Shore J. Stat. Phys. 78, 1277 (1995). 



12 




Figure 1: The six vertex configurations and their related Boltzmann weights. 
In (a) we draw all the arrows and in (b) we draw by solid lines th links where 
the arrows are pointing to the down and left directions. 
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